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Abstract 

We consider the hexagonal circle packing with radius 1/2 and perturb it by letting the circles 
move as independent Brownian motions for time t. It is shown that, for large enough t, if Il t 
is the point process given by the center of the circles at time t, then, as t — > oo, the critical 
radius for circles centered at II t to contain an infinite component converges to that of continuum 
percolation (which was shown — based on a Monte Carlo estimate — by Balister, Bollobas and 
Walters to be strictly bigger than 1/2). On the other hand, for small enough t, we show (using 
a Monte Carlo estimate for a fixed but high dimensional integral) that the union of the circles 
contains an infinite connected component. We discuss some extensions and open problems. 

1 Introduction 

Let T be the triangular lattice with edge length 1 and let IIo be the set of vertices of T ■ We see 
no as a point process and, to avoid ambiguity, we use the term node to refer to the points of no- 
Now, for each node u £ Hq, we add a ball of radius 1/2 centered at u, and set R(Jlo) to be the 
region of M? obtained by the union of these balls; more formally, 

R(U ) = (J B(x,l/2), 
xerio 

where B(y,r) denotes the closed ball of radius r centered at y. In this way, -R(no) is the so-called 
hexagonal circle packing of R 2 ; refer to [I] for more information on packings. Clearly, the region 
R(Ho) is a connected subset of R 2 . 

Our goal is to analyze how this set evolves as we let the nodes of Hq move on R 2 according 
to independent Brownian motions. For any t > 0, let H be the point process obtained after the 
nodes have moved for time t. More formally, for each node u € no, let (Cu(*))t be an independent 
Brownian motion on R 2 starting at the origin, and set 

n* = |J (u + ( u (t)). 

iterio 

A natural question is whether there exists a phase transition on t such that R(Ilt) has an infinite 
component for small t but has only finite components for large t. 

Intuitively, for sufficiently large time, one expects that will look like a Poisson point process 
with intensity 2/y/3, which is the density of nodes in the triangular lattice. Then, for sufficiently 
large t, R(Jlt) will contain an infinite component almost surely only if R(&) contains an infinite 
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component where $ is a Poisson point process with intensity A = 2/\/3. In the literature, R(&) 
is referred to as the Boolean model. For this model, it is known that there exists a value A c so 
that, if A < A c , then all connected components of R(&) are finite almost surely [9]. On the other 
hand, if A > A c , then i?(3>) contains an infinite connected component. The value of A c is currently 
unknown and depends on the radius of the balls in the definition of the region R. When the balls 
have radius 1/2, it is known that A c satisfies 0.52 < A c < 3.38 [3] Chapter 8], but these bounds 
do not answer whether 2/\/3 is smaller or larger than A c . However, using a Monte Carlo analysis, 
Balister, Bollobas and Walters [2] showed that, with 99.99% confidence, A c lies between 1.434 and 
1.438, which are both larger than 2/y/3. We then have the following theorem, whose proof we give 
in Section [3] 

Theorem 1.1. If A c > 2/y/3, where A c is the critical intensity for percolation of the Boolean model 
with balls of radius 1/2, then there exists a positive to so that, for all t > to, R(Ht) contains no 
infinite component almost surely. 

We now give a different statement of Theorem ll.il For each t, there exists a critical radius r c (t) 
so that, adding balls of radius r > r c (t) centered at the points of IL gives that the union of these 
balls contains an infinite component almost surely. Then, the next theorem, whose proof we give 
in Section [3] follows from the proof of Theorem 11.11 and a result of Sinclair and Stauffer [13] . 

Theorem 1.2. As t — > oo, we have that r c (t) converges to the critical radius for percolation of the 
Boolean model with intensity 2/\f?>. 

In order to prove Theorems 11.11 and 11.21 we devote most of Section [3] to prove Theorem 11.31 
below, which we believe to be of independent interest. Consider a tessellation of R 2 into regular 
hexagons of side length 5y/i, where 5 is an arbitrarily small constant. Let I denote the set of points 
of R 2 that are the centers of these hexagons. Then, for each i E I, denote the hexagon with center 
at i by Qi and define a Bernoulli random variable AQ with parameter p independently of the other 
Xj, j i. Define 

c( P ,s)= U Qi- 

iel: X i= l 

When p > 1/2, which is the critical probability for site percolation on the triangular lattice, we 
have that C(p,5) contains a unique infinite connected component. We are now ready to state our 
main technical result, whose proof is given in Section [3] 

Theorem 1.3. There exists a universal constant c > and, for any p £ (0, 1) that can be arbitrarily 
close to 1 and any arbitrarily small 5 > 0, there exists a to > such that, for all t > to, we can 
couple lit, (Xi)i£i and a Poisson point process of intensity + c\fb, which is independent of 
Wo, so that 

R(u t )nc(p,5) c R($)nC(p,S). 

In words, Theorem [L3] establishes that lit is stochastically dominated by a Poisson point process 
inside C(p, 5). As we show later in Lemma [2.1 1 LTj cannot be stochastically dominated by a Poisson 
point process in the whole of M 2 . 

We note that the opposite direction of Theorem 11.31 was established by Sinclair and Stauffer |13[ 
Proposition 4.1], who proved that, under some conditions on the initial location of the nodes, after 
moving as independent Brownian motions for time t, the nodes stochastically dominate a Poisson 
point process. The result of Sinclair and Stauffer has been used and refined in [12], [TJ], and turned 
out to be a useful tool in the analysis of increasing events for models of mobile nodes, such as 
the so-called percolation time in [12j and detection time in [T4]. We expect that the ideas in our 
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proof of Theorem 11.31 can help in the analysis of decreasing events, which have so far received less 
attention. 

We remark that the proof of our Theorem 11.31 requires a more delicate analysis than that of 
Sinclair and Stauffer. In their case, nodes that ended up moving atypically far away during the 
time [0, t] could be simply disregarded as it is possible to show that the typical nodes already 
stochastically dominate a Poisson point process. In our setting, no node can be disregarded, 
regardless of how atypical its motion turns out to be. In order to solve this problem, we first 
consider what we call well-behaved nodes, which among other things satisfy that their motion 
during [0, t] is contained in some ball of radius c\/t for some large constant c (we defer the complete 
definition of well-behaved nodes to Section [3|). The definition of well-behaved nodes is carefully 
specified so that any given node is likely to be well-behaved and, in addition, it is possible to show 
that well-behaved nodes are stochastically dominated by a Poisson point process inside C(p, S). For 
the remaining nodes, which comprise only a small density of nodes, we use a sprinkling argument 
to replace them already at time by a (non-homogeneous) Poisson point process of low intensity. 
Then, even though the motion of the nodes that are not well behaved is hard to control, we use 
the fact that they are a Poisson point process at time to show that, at time t, they stochastically 
dominate a Poisson point process of low intensity. 

Now, for the case when the nodes of Ho move for only a small time t, we believe the following 
is true. 

Conjecture 1.4. There exists a to > so that, for all t < to, -R(IL) contains an infinite component 
almost surely. 

We are able to establish the conjecture above given a Monte Carlo estimate for a finite but high 
dimensional integral. We discuss the details in Section HI Note that if Conjecture 11.41 is true, then 
a curious consequence of this and Theorem 11.21 is that r c (t) is not monotone in t. 

We conclude in Section [5] with some extensions and open problems. 



2 Stochastic Domination 

We devote this section to the proof of our main technical result, Theorem 11.31 where we study the 
behavior of the balls after they have moved for a time t that is sufficiently large. In Section [3] we 
show how to adapt this proof to establish Theorems 11.11 and 11.21 

Intuitively, as t — > oo, 11^ looks like an independent Poisson point process of intensity 2/\/3. 
Since the intensity of $ is larger than 2/\/3, we would like to argue that there exists a coupling 
between <3? and IT such that contains IT. Unfortunately, this cannot be achieved in the whole 
of M 2 , as established by the lemma below, which gives that, for any fixed t, the probability that a 
sufficiently large region 5cl 2 contains no node of lb is smaller than exp ^— ^-^= + cy/d^j vol (S)j , 
which is the probability that <E> has no node in S. Hence, cannot stochastically dominate lit m 
the whole of R 2 . 

Lemma 2.1. Fix t sufficiently large and let S be a hexagon of side length (log t)\/t obtained as the 
union of (log t) 2 t triangles ofT- Then, there exists a positive constant d such that 

P (lit n S = 0) < exp (-c'(log t) 2 vol (S)) . 

Proof. For simplicity we assume that (log t)y/i is an even integer. Let x be the middle point of S 
and consider the hexagon S' of side length ^^y/i composed of ( lo |^ t triangles of T and centered 
at x. Note that S contains the ball B (x, v ^; og * \/t ) and 5' is contained in the ball B (x, ^y^-y/i ) • 
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Therefore, a node of IIo that is inside S' can only be outside of S at time t if it moves at least 
— — ^— —y/i > ^p-y/i. For any fixed node u G rig, we have from the Gaussian tail bound (cf. 
Lemma |A.3|) that 

logt /-\ 3 / {\ogtf 



18 



Each node of LTo belongs to 6 triangles of T, then there are at least ^°J^ t nodes of LIo in S' . Since 
each of them need to move more than ^&y/t by time t for S to contain no node of Il t , we obtain 



P(n t n5 = 0)< -= exp - v ° ' < exp 



(logt) 2 \\ « * / (log*) 4 * 



^Trlogt V 18 /y " V 432 

Since vol (S) = ^p(logi) 2 i, the proof is completed. □ 

Now we turn to the proof of Theorem 11.31 The goal is to show that 3> contains IT inside a 
percolating cluster of a suitable tessellation of R 2 . For this, we tessellate IR 2 into hexagons of side 
length 8\fi. We take this tessellation in such a way that no point of LTo lies on the edges of the 
hexagons; this is not crucial for the proof but simplifies the explanations in the sequel. Let H. 
denote the set of these hexagons. Consider a node v £ LTo. Let Qi be the hexagon of H. that 
contains v and let v' be a copy of v located at the same position as v at time 0. We let v' move 
up to time t according to a certain procedure that we will describe in a moment, and then we say 
that v is well behaved if we are able to couple the motion of v with the motion of v' so that v and 
v 1 are at the same location at time t. Recall that I is the set of points given by the centers of the 
hexagons in T~L. For i E I, we define 



J, 

where 



j € I: sup ||x — y\\2 < C8\/i >, 



c = 4r 3 / 2 . 



For i,j such that j E Jj we say that i and j are neighbors. Now we describe the motion of v' . Let 
ft be the density function for the location of a Brownian motion at time t given that it starts at 
the origin of ]R 2 . We fix t and, for each i,j £ I such that i and j are neighbors, we let 

<Pt(hj) = i nf „ ft(y-x). (1) 

If i and j are not neighbors we set ipt(i,j) = 0. Then, the motion of v' is described by first choosing 
a j G Jj with probability proportional to (ft(i,j) and then placing v' uniformly at random in Qj. 
The main intuition behind this definition is that, when v is well behaved, its position inside Qj has 
the same distribution as that of a node of a Poisson point process inside Qj. Therefore, as long 
as the number of well behaved nodes that end up in Qj is smaller than the number of nodes in 
$ n Qj, we will be able to couple them with <£. Another important feature for the definition of well 
behaved nodes is that, if v is well behaved and ends up moving to hexagon Qj, then we know that, 
at time 0, v was in some hexagon of Jj. In particular, there is a bounded number of hexagons from 
which v could have moved to Qj, which allows us to control dependences. 

Now we show that nodes are likely to be well behaved. Since the area of each hexagon of H is 
'^■5 2 t, we have that 

3\/3 

P (y is well behaved) = - 5 2 t(p t (i,j). (2) 

jeJi 
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The idea is that 5 is sufficiently small so that ft varies very little (i.e., ft is essentially constant) 
inside any given hexagon of H, but, at the same time, C5 is large so that the probability that v 
moves to an hexagon that is not in J{ is small. We can then obtain in the lemma below that the 
probability that v is well behaved is large. 

Lemma 2.2. Let v be a node ofHo located in Qi. We have 

(C-3) 2 < |^| < |c 2 , 

and, for sufficiently large t, we have 

P (v is well behaved) > 1 — 55. 

Proof. For j G" Jj, we know, by definition, that there exist a xq G Qi and a y§ G Qj such that 
ll^o — J/o II 2 > CSy/i. Then, by the triangle inequality, we have that, for any y G Qj, 

\\y - * || 2 > Cb\ft - \\y - 2/0 1| 2 - IN ~ x oh > CSVt - 35Vt, (3) 

where we used the fact that, for any two points y, yo in the same hexagon, we have \\y — j/o II2 < 2b\ft 
and, for any x G Qi we have \\i — 1 1 2 < o~Vi. Therefore, if we add balls of radius 5\/t centered at 
each j G Jj, these balls cover the whole of B{i, CSy/i — 35^), which yields 

vol{B(i,C5Vt-35Vt)) _ o ^ 2 
W - vol(B(0,M)) " ( J ' 

For the other direction, note that if we add balls of radius ^Sy/i centered at each j G J%, these 
balls are disjoint and their union is contained in B{i,Cb\ft), which gives 

vol {B(z,C8Vt)) 4 2 
vol f B (0, ifSVi 

Now we prove the second part of the lemma. Note that, using ([2]) and ([I]), we have 
P (v is well behaved) = y, - — exp 

>~ E/ Q ^e Xp (-fc^W) rf2 , (4) 
jeJi JQ j v 7 

where the last step follows by the triangle inequality. Now, from Q, the ball B(i,CS\/t — 35^/t) 
only intersects hexagons that are neighbors of i. We denote by S a = [—a/2, a/2] 2 the square of 
side length a, and, for any z = (z\,Z2) G M 2 and a G R+, we use the inequality (||z||2 + a) 2 < 
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(|zi| + a) 2 + (|z2| + a) 2 ■ Then, applying Q, we obtain 
P (v is well behaved) > j - — exp 



> 



> 



> 




5 2 (C + 3(V2-1)) 2 \\ 2 



where the second to last step follows by the standard Gaussian tail bound (cf. Lemma lA.3p . Then, 
using the value of C, we have that ^ g|? e~ <52c ' 2//4 < -^e -4 for all 5 £ [0, 1]. Then, the right-hand 

side above is at least 1 — 5 ( -7= + ^7= ) > 1 — 55. □ 

V V2?r V 71 " / ~~ 



We will treat the nodes that are not well behaved by means of another point process. For any 
point x 6 IR 2 , we set q(x) = i if x € Qi- Then, let gt{x, y) be the density function for a node v that 
is not well behaved to move from x to y after time t. We have that 

/ \ ft(y ~ x) - ip t (q(x),q(y)) 

9t(x,y) = WT~- I rTTT Tr ( 5 ) 

P [v is not well behaved) 

For each v £ IIq, let N v (fi) be a Poisson random variable with mean /i, and let ^o(fJ-) be the point 
process obtained by putting N v (fi) points at v for each v £ ILj. We set e _At = P is well behaved) 
and, from Lemma [2 .21 and the fact that 5 is sufficiently small, we henceforth assume that (i < 1. We 
can then use a standard coupling argument so that N v (/j,) > 1 if and only if v is not well behaved. 
The intuition is that, by replacing each node of LTo that is not well behaved by a Poisson number 
of nodes, we can exploit the thinning property of Poisson random variables to show that, as the 
nodes move, they are stochastically dominated by a Poisson point process. 

For each w £ ^oO-O; let £ w (t) be the position of w at time t according to the density function 
g t . Define ^t(^) to be the point process obtained by 

*t(M)= (J &»(*)• 

The following lemma gives that VPf(/i) is stochastically dominated by a Poisson point process. 

Lemma 2.3. Let be the probability that a node of Uq is well behaved. For t sufficiently large, 
there exists a universal constant c > such that, i/ \E r is a Poisson point process with intensity c\fb, 
then it is possible to couple with ^(/i) so that C ty. 

Proof. Since the nodes of &o(n) move independently of one another, we can apply the thinning 
property of Poisson random variables to obtain that ^t(n) is a Poisson point process. Let A(x) 
be the intensity of ^t(lJ-) at x £ M?. By symmetry of Brownian motion and the symmetry in the 
motion of well behaved nodes, we have that 

A ( x ) = ^9t(v,x) = ^ m(x,v). (6) 
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Recall that, for any z£l 2 and £ > 0, we define z + Se as the translation of the square [0, £] 2 so 
that its center is at z. Define the square R\ as x + S 5S ^q, the annulus R 2 as [x + S 5CS ^q) \ R\ and 
the region R% as M 2 \ (Ri U R2). We split the sum in ((6|) into three parts by considering the set of 

points Pi = n n Ri, p 2 = n n r 2 and p 3 = n n i? 3 . 

We start with P 2 . We can partition each hexagon of T~L into smaller hexagons of side length 
y/3/3 such that each point of IIo is contained in exactly one such hexagon. This is possible since the 
dual lattic43 of T is a hexagonal lattice of side length y/3/3, so the hexagons of side length y/3/3 
mentioned above can be obtained by translating and rotating the dual lattice of T ■ We denote by 
H 1 the set of hexagons of side length y/3/3 obtained in this way. 

For each z € M 2 , let H z be the hexagon that contains z in T~L' . Each H z has side length 
V3/3 and area y/%/2. Thus, for any point z G R 2 , we have that H z C x + S$c5-Ji+W3/Z' Let 
#2 = (x + S 5C8V - t+4V z /3 ) \ Ri, which gives that 

^ Hg t (x,v) < —= / sup fig t (x,z')dz. 



veP 2 v " ^2 



V3 



ill, z'e-ffz 



Now, note that p(t; is not ^ cll behavcd) = < < 2 since p < 1. Then, using the definition 

of from ([5]) and the definition of ipt in ([!]), we have that 



fig t (x,v) < [ ( sup f t (z - x) - ip t {q{x),q{z')))dz 

~B V3 Jr> z'eHz 



veP 2 v " ^2 



- ~at / ( SU P M z> -x)- inf /((a/ - z") ) cfe. 

V3 7^ \z'eH z x'£Q q ( X ),z"£Q q{z >) 



Now, by the triangle inequality, we have that \\z' — x\\ 2 > \\z — x\\ 2 — \\z — z'\\ 2 > \\z — x\\ 2 — 2y3/3 
and \\x' — z"\\ 2 < \\z — x\\ 2 + \\z — z"\\ 2 + ||x — x'\\ 2 < \\z — x\\ 2 + A5y/i, which gives that 

v£P 2 



Note that we can write 



{\z-X2- 
exp : — — l — exp 




| 2 + W*) : 



2t 



2 sVi + 88 2 ty 



exp — exp exp 



t \ It 



Now we use that, for z € R' 2 , we have \\z — x\\ 2 < 5v ^ g Sy/i + 2\/6/3. Then, the first exponential 
term above is 1 + o(l) and, for the second exponential term, we can use the inequality e~ x > 1 — x, 
which gives, as t — > 00, 

«) < ^( 10 ^ 2 + 852 + o ( 1 ))/ J ,^ ex p(-^ir^) dz 



v£P 2 



< Cl VS + o{l), (7) 



1 Recall that the dual lattice of T is the lattice whose points are the faces of T and two points are adjacent if their 
corresponding faces in T have a common edge. 
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for some universal constant c\ > 0. 

For the terms of © where v G P3 we have that gt(x,v) = p ^ is not* wolf behaved) ' Then, let 
R' 3 = M 2 \ (x + SscsVi-4\/3/3^ so that, for each z G i?3, we have il 2 C -R' 3 , which allows us to write 



V I \s 4 I II >\J s 4 f 1 ^ (lk-»ll2-2^/3) 2 \ 

4 s 75 i, ss. /t(I - 2 » rf2 - vs 4 2s exp ( » J dz - 

Now, letting to = z — x and writing w = (w\,W2) we have that 

(\\w\\ 2 - 2v / 3/3) 2 > (K| - 2v / 3/3) 2 + (> 2 | - 2^/3)2 - 4/3, 
which can be used to get the bound 

V ugt(x,v) < 4=exp (— M 2 / exp dtoi | <^r+o(l), (8) 

for some universal constant C2 > 0. Finally, for the terms in Q with v G Pi, we use that 
figt(v,x) < 2ft(v,x) < for all v,x which gives that 

^ m ( x , v )<±^(^5Vt + 4^pJ <c 3 5 2 + o(l), (9) 

for some universal constant C3 > and where ^56y/i + 4-^^ is an upper bound for the number 
of points in P\. 

Plugging ©, © and © into © yields 

A(z) = m{x,v) < 01^+^ + 03^ + 0(1). 

□ 

We now proceed to the proof of Theorem 11.31 

Proof of Theorem 11.31 We start by giving a high-level overview of the proof. First, we assume 
that all nodes of Ilo are well behaved. Then we consider a hexagon Qi of T~L and count the number 
of such well behaved nodes that are inside Qi at time t. Note that, by the definition of well behaved 
nodes, given that a node is in Qi at time t, then its location is uniformly random in Qi. Therefore, 
in order to show that they are stochastically dominated by a Poisson point process, it suffices to 
show that there are at most as many nodes of Ilo in Qi at time t as nodes of the Poisson point 
process. This will happen with a probability that can be made arbitrarily large by setting t large 
enough. We then use the fact that, since nodes are considered well behaved, a node can only 
be in Qi at time t if that node was inside a hexagon of Jj at time 0. Therefore, if we consider 
a hexagon Qj such that Jj n Jj = 0, we have that the well-behaved nodes that are able to be 
in Qi at time t cannot end up in Qj. Hence, the event that the well-behaved nodes in Qi are 
stochastically dominated by a Poisson point process is independent of the event that the nodes in 
Qj are stochastically dominated by a Poisson point process. This bounded dependency is enough 
to complete the analysis of well behaved nodes. On the other hand, to handle nodes that are not 
well behaved, we add a discrete Poisson point process at each node of Ilo so that the probability 
that we add at least one node at a given u G Ilo is exactly the same as the probability that v is 
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not well behaved. Thus, this discrete Poisson point process contains the set of nodes that are not 
well behaved. We then use Lemma 12.31 to conclude that the nodes that are not well behaved are 
stochastically dominated by a Poisson point process, which concludes the proof. 

We now proceed to the rigorous argument. For each v £ Ho, let £' v (t) be the position of v at 
time t given that v is well behaved, and let 

n; = |J m- 

ueiio 

Note that, since e _M is the probability that a node is well behaved and ^o(^) 1S the point process 
obtained by adding a random number of nodes to the points of IIo according to a Poisson random 
variable with mean /i, then there exists a coupling so that 

n f cn|uf t (4 

Lemma 12.31 establishes that ^(/x) is stochastically dominated by a Poisson point process with 
intensity c\y§ for some universal constant c\ > 0. It remains to show that Ii' t is also stochastically 
dominated by a Poisson point process. Unfortunately, this is not true in the whole of M? as shown in 
Lemma l2.1i We will then consider the tessellation given by H and show that, for each hexagon Qi 
of the tessellation with Xi = 1, where the Xi are defined in the paragraph preceding Theorem 11.31 
U' t is stochastically dominated by a Poisson point process II of intensity (1 + y/5)2/V3. 

In order to see this, for each i £ /, we define a binary random variable Yj, which is 1 if LT has 
more nodes in Qi than Tl' t . Then, since each node of U' t is well behaved, whenever Yj = 1, we can 
couple LT with H' t such that IT 2 n£ in Qi. First we derive a bound for the number of nodes of H' t 
inside Qi. For each v £ IIo, let Z v be the indicator random variable for £' v (t) £ Qi. Then, since the 
probability that £' v (t) £ Qi is proportional to ipt(q(v),i), the expected number of nodes of W t in Qi 
is 

y e[z v }= y ^M= y f^R = 3S % 

^ L J ^ M ^ M 

veiLoniu^j.Qj) t>6non(u ieJi Qj) t)enon(Uj eJi Q^) 

where M is a normalizing constant so that ^ - <ft(i,j) = M for all i. The last step follows since 
Su6n n(Uj 6 j.Qj) Pt(i,q{v)) = 35 2 t X^jeJ, Ptihj) = 35 2 tM. A simpler way to establish the equation 
above is by using stationarity and noting that 35 2 t is the number of points of IIo i n Qi- Since 
the random variables Z v are mutually independent, we can apply a Chernoff bound for Binomial 
random variables (cf. Lemma |A.2|) to get 

\^en n(u ieJ .Q,-) / \ 



where the last step follows from Lemma 12.21 Using a standard Chernoff bound for Poisson random 
variables (cf. Lemma I A. 1|) we have 

P (il has less than (1 + \f5/2)?>5 2 t nodes in Qi) < exp 



2(1 + ^)7 ' 



Therefore, we obtain a constant C2 such that 



P(Y = l)>l-exp(-^). (10) 

The random variables Y are not mutually independent. However, note that Yi depends only on the 
random variables Yy for which Jj/ n Jj ^ 0. This is because, for any i £ I, only the nodes that are 
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inside hexagons Qj with j £ Jj can contribute to Yj. Therefore, using Lemma 12.21 we have that Yj 
depends on at most (|C 2 ) 2 other random variables Y . By having t large enough, we can make the 
bound in (jlOp be arbitrarily close to 1. This allows us to apply a result of Liggett, Schonmann and 
Stacey [7J Theorem 1.3], which gives that the random field (Yj)j e j stochastically dominates a field 
(Y-)i £ j of independent Bernoulli random variables satisfying 

p(y/ = i)>i-ex P (-^r), 

for some positive constant C3. So, with t sufficiently large, we can assure that P (Y- = 1) is larger 
than p in the statement of Theorem 11.31 Then, we have that, whenever Y- = 1, the Poisson point 
process nuvt'tC/-*) stochastically dominates IT inside Qj. Since II and ^t(^) are independent Poisson 
point processes, we have that their union is also a Poisson point process of intensity no larger than 

-l + ^ + cA (11) 

which completes the proof of Theorem 11.31 □ 



3 Large time 

In this section we give the proofs of Theorems 11.11 and 11.21 Here we will use some steps of the proof 
of Theorem 11.31 without repeating the details. For this reason, we suggest the reader to read the 
proof of Theorem 11.31 before embarking in the proofs in this section. 

Proof of Theorem 11.11 For each i £ I, let TV, be the set of hexagons Qj such that Qi and Qj 
intersect. Now, define a binary random variable Yi to be 1 if Yj = 1 for all j £ Aj and the largest 
component of R((Uj£]y i Qj)Pi (LTU^t (//))) has diameter smaller than 5y/t/W. (Recall the definition 
of Yi from the paragraph preceding (jlOp .) If A c > we can set 5 small enough so that the 

intensity of n U ^(/z), which is bounded above by (llip . is smaller than A c . Then, using (jlOD and a 
result of Penrose and Pisztora [111 Theorem 2], we have that, for any given i and sufficiently large 
t, 

P (Yi = l) > 1 - 6exp " exp (-c 4 <Vi) , (12) 

for some universal constant C4 > 0. Also, the Yj depends only on the random variables Yy for 
which Ji' n Jj 7^ for all j £ iVj. Hence, Yj depends on no more than 7(|C 2 ) other Yj since 
\Ni\ = 7. Thus, as in the proof of Theorem 11.31 we can a PPb/ the result of Liggett, Schonmann 
and Stacey Theorem 1.3] to conclude that (Yi)i stochastically dominates a random field of 
independent Bernoulli random variables with mean p, which can be made arbitrarily close to 1 
by having t large enough. As a consequence, we have that, almost surely, the set of hexagons 
with Yi = has only finite components. With this, we can conclude that R(Ht) has no infinite 
component almost surely since the coupling between R^, (Yi)j, n and ^ft(^) gives that a connected 
component of R(Ht) can only intersect two non-adjacent edges of an hexagon i £ / if Yj = 0. 
But we showed that any set of intersecting hexagons with Yj = is finite and, therefore, must be 
surrounded by hexagons j with Yj = 1, which are not crossed by R(Jlt). Hence, all components of 
R(Ht) are finite. □ 

Proof of Theorem 11.21 In this proof we will use some steps from the proof of Theorem 11.11 Let 
be the critical radius for percolation of the Boolean model with intensity A. Therefore, if t is 
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large enough and we add balls of radius r < r^° centered at the points of IT, where Ao is the 
intensity of IT U ^t(/i), which is bounded above by (jlip . we have that (|12p holds and the union of 
the balls do not have an infinite component. Now, we have that the Boolean model with intensity 
A and radius r is equivalent (up to scaling) to the Boolean model with intensity A and f provided 
Ar 2 = Af 2 . Therefore, for any e > 0, we have that 



J^r e (13) 



Using this we obtain that, for any 5 > 0, 



lim inf r c (t) > r c = r 



A„ _ r 2/v / 3 i 



where c > is a universal constant. Since 5 can be arbitrarily close to 0, we obtain 

liminf r c (i) > r 2/v ^. (14) 



t— >oo 



Now, to obtain an upper bound for r c (t), we use a result of Sinclair and Stauffer [13|. Proposi- 
tion 4.1], which can be stated as follows. 

Proposition 3.1 ( \13\ Proposition 4.1]). Consider a square S of side length K tessellated into 

subsquares of side length I, and assume that each subsquare contains at least (3£ 2 nodes at time 

0. Denote the nodes by Ho and let Ha be the point process obtained after the nodes of Ho have 

moved as independent Brownian motions for time A . Then, there are positive constants c\ , C2 , C3 

such that, for any e > 0, if A > c\l 2 /e 2 , we can couple Ho with an independent Poisson point 

process of intensity (1 — e)/3 so that, inside a square S' of side length K' < K — C2 \/A log e _1 with 

the same center as S, Ha contains the nodes of the Poisson point process with probability at least 
1 _ e -c 3 e 2 /3e 2 _ 

Now we show how to apply this result to our setting. Let i £ I be fixed. Take S to be the union 
of the hexagons in J^; clearly S is not a square but that is not important. Instead of tessellating S 
into squares, we tessellate S into hexagons of side length \/3/3. Each such hexagon contains one 
node of Ho = Ilo n (Uj e j a (5j), which gives f3 = 2/\/3, the density of nodes of Ho- The main step 
in adapting the proof to our setting is to note that I can represent the diameter of the cells of the 
tessellation, which in our case gives I = 2\/3/3, and K - K' is the minimum distance between S' 
and a point outside S. Since we take S' = Qi we have that K — K' > C5y/t - 25Vt, where 25^/t 
is the diameter of Qi. Now setting e = Vd and A = t, we obtain a to so that, for all t > to, the 
conditions on A and K — K' in Proposition 13.11 are satisfied. Therefore, with probability at least 
1 — e _C3<53 *, for some positive constant C3, the nodes of IT that are inside Qi at time t and were 
inside Uj 6 j 4 Qy at time stochastically dominate a Poisson point process $ of intensity (1 — ^/5)-^. 
When this happens, we let Yi = 1, where Yi here is analogous to the one in (llOp . 

Then we proceed similarly as in the proof of Theorem II. li Define Yi to be 1 if Yj = 1 for all 
j € Ni (recall that iV, is the set of hexagons that intersect Qi) and the largest component of the 
region R(^)n(Uj^NiQj), which we denote by X, is such that Uj£NiQj\X contains only components 
of diameter smaller than 5\/t/10. This means that if there exists a path Y\,Y2, . . . such that, for all 
j = 1, 2, . . . , we have Yj = 1 and Qj and Qj+i intersect, then the region R(&) n (Q\ U Q% U ■ ■ ■ ) has 
a connected component that intersects each hexagon Qj in the path. Then, for t sufficiently large, 
we can apply [111 Theorem 2] and Theorem 1.3] as before to show that the hexagons with Yi = 1 
percolate and, consequently, adding balls of radius r > r c centered at the nodes of $ produces 
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an infinite component. By the scaling argument in (|13p we have that S = T\J^ \ T\ 

which finally yields 



hmsupr c (t) < r 2 J^-^ = ^J—*—. 

t-^oo y A — v 5 

Since 5 can be arbitrarily close to 0, we obtain lim sup^^ r c (t) < r% 3 , which together with (fl4 
concludes the proof of Theorem 11.21 □ 



4 Short time 



Now we turn our attention to the case when t is sufficiently small. We establish that, given a Monte 
Carlo estimate, R(Ht) contains an infinite component almost surely for sufficiently small t. 

Consider a tessellation of M. 2 into regular hexagons of side length 50. We will denote this 
tessellation by Hso- Instead of considering the usual tessellation, where each hexagon is obtained 
by the union of some triangles of T, we will shift the hexagonal tessellation (see the illustration 
in Figure [T]) so that no node of n is on an edge or vertex of H50, and the edges of % 50 intersect 
as many of the balls centered at Hq as possible. More formally, since a transitive lattice can be 
specified by a single edge, we define T as the triangular lattice containing an edge between the 
points (0,0) and (1,0), and for any I > 0, we let Tin be the hexagonal lattice containing an edge 
between (1/2, -\/3/4) and (£ + 1/2, -y/3/4). 

Let Hi and H2 be two hexagons of T~L^q that have one edge in common, and denote this edge 
by e. Starting from e, denote the other edges of H\ in clockwise direction by e%, &21 63, e^, es; thus 
e% is the edge of H\ opposite to e. Similarly, denote the other edges of H2 in clockwise direction 
by e' 1; e' 2 , eg, e^, e' 5 (refer to Figure [I]). Given any three sets X\, X 2 , X3 C M 2 and any t > 0, we say 
that R(TLt) has a path from X\ to X2 inside X$ if there exists a sequence of nodes ui,«2, ■■■ ,Uj- 
of lb, all of which inside X%, such that B{u\, 1/2) intersects X±, B(uk, 1/2) intersects X2, and for 
each i > 1, the distance between m and Uj+i is at most 1. With this, we say that R(Ht) crosses 
Hi and H2 if the following three conditions hold: 

1. R(Ht) has a path from 63 to e 3 inside Hi U i?2- 

2. R(Ht) has a path from e\ U e2 to U es inside iTi U i?2- 

3. -R(ILJ has a path from U e 2 to e^ U e^ inside Hi U i?2- 

We denote by At the event that R(Ht) crosses Hi and H2 with a path that also crosses Hi and 
i?2 at time 0; this last condition is used to obtain a type of monotonicity later. Then we have the 
following result. 

Theorem 4.1. Suppose that there exists an eo > such that P (A t0 ) > 0.8639. Then, for all 
e G [0, eo], ZioWs that -R(n e ) contains an infinite connected component almost surely. 

Note that, for any fixed t, verifying the condition P (At) > 0.8639 resorts to solving a finite, but 
high dimensional integral describing the crossing probability. We were able to check the validity of 
this condition for t = 0.01 via a Monte Carlo analyst with confidence 99.99%. 

We start the proof with the lemma below, which establishes a type of monotonicity that will be 
useful later. To state the lemma, let V be a set of points in M rf and define E(V) as the set of pairs 
of points of V whose distance is at most 1. Note that the pair (V, E(V)) induces a graph over V. 

2 To obtain this Monte Carlo estimate we employed the Mersenne Twister pseudorandom number generator by 
Matsumoto and Nishimura [8] with period 2 19937 — 1 and improved initialization scheme from January 26 th , 2002. 



12 





J>— •e 4 -' : " /" 







































.Jft A-e- 









Figure 1: The hexagonal tessellation H4 of M? with hexagons of side length 4, and the illustration 
of a path of intersecting balls crossing Hi and H2 ■ 



Lemma 4.2. Let V C IIo sitc/i £/iai i/ie graph (V,E(V)) is connected. Let V s be obtained by letting 
the nodes of V move for time s > according to independent Brownian motions. Then, we have 
that P (E(V) C i?(V^)) is non-increasing with s. 

Proof. This follows by Brownian scale. Consider s' > s, and let V', be obtained by letting the 
nodes of V move according to independent Brownian motions for time s' . Now we create a coupling 
between £(V S ',) and E(V S ) so that, if E(V) C S(V^), then C 

For each let C«( s ) an d C'u( s ') be the Brownian motions for the motion of u in V s and VI'/, 

respectively. Then, 

P(E(V)QE(V s )) = pI fl {||u + Cu(8)-t>-C«(*)l|2<l} J. 

\(«,«)es(y) / 

Now, by Brownian scale, we can couple Cu(s) and C' u ( s ') v ^ a Cu(s) = \/ $/ s'C u ( s ')- Using this, we 
write the right-hand side above as 

n {\\ u -v+^((' u (s')-as'))h<i}) . 

Now define the vectors x\ = u — v and %i = C u (s') — Cv( s ')- For any 5 E (0, 1), it follows by standard 
geometric arguments that, if 1 1 acx 1 1 2 an d ||xi + 1| 2 are at most 1, then 



Ixi + 6X2W2 < 1, 



which establishes that 
P 



> P 



{\\u-v+ yfip{&M ~ as'))h < l} I 
u,v)eE(V) J 



Ju,v)eE(V) 



>) 
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which completes the proof. 



□ 



Now we proceed to the proof of Theorem 14.11 

Proof of Theorem 14.11 Note that, due to Lemma 14.21 proving Theorem 14.11 reduces to showing 
that, if there exists an e > such that P (A e ) > 0.8639, then R(H e ) contains an infinite connected 
component almost surely. We henceforth fix a value of e and assume that P (A e ) > 0.8639. 

We will use a renormalization argument. Consider the hexagons T-L^q described in the beginning 
of this section. Now, define the graph L = (U, F) such that U is the set of points given by the 
centers of the hexagons and F is the set of edges between every pair of points i,j G U for which the 
hexagons with centers at i and j share an edge. Note that L consists of a scaling of the triangular 
lattice. 

We now define a collection of random variables Xi for each edge i € F. In order to explain the 
process defining Xi, let H\ and H2 be the hexagons whose centers are the endpoints of i. We then 
define Xi = 1 if and only if .R(n e ) crosses H\ and H2 with a path of balls that also crosses H\ and 
H2 at time 0. (The definition of crossings is given right before the statement of Theorem 14.11 ) Let 
j be an edge such that % and j are disjoint, and let #3 and H± denote the hexagons centered at 
the endpoints of j. Clearly, Xi and Xj are independent since the set of balls crossing H\ and H2 
at time does not intersect the set of balls crossing H3 and at time 0. Thus, the collection 
(Xi)i is a so-called 1-dependent bond percolation process, with P (Xi = 1) = P (A t ) > 0.8639. 
Then, we can use a result of Balister, Bollobas and Walters [H Theorem 2], which gives that any 
1-dependent bond percolation process on the square lattice with marginal probability larger than 
0.8639 percolates almost surely. Since the triangular lattice contains the square lattice, we obtain 
that, almost surely, there exists an infinite path of consecutive edges of F with Xi = 1 for all i in 
the path. 

To conclude the proof, note that, for two non-disjoint edges i and j with Xi,Xj = 1, we have 
that the crossings of the hexagons whose centers are located at the endpoints of i and j intersect. 
Then, the infinite path of Xi with Xi = 1 for all i in the path contains an infinite path inside 
R(IL e ), which concludes the proof of Theorem 14.11 □ 

5 Extensions and Open Problems 

In the remaining of this section we discuss extensions and open problems regarding other circle 
packings (Section 15. ip . balls moving over graphs (Section \5.2\i and critical radius for non- mobile 
point processes (Section I5.3P 

5.1 Other circle packings 

Let LTq be the point process given by the vertices of the square lattice with side length 1, and let 
be the point process obtained by letting the nodes of LTq move for time t according to independent 
Brownian motions. Note that, for any e > 0, if we look at two balls of radius 1/2 centered at 
two adjacent nodes of LTq, then at time e, the probability that these two balls intersect is strictly 
smaller than 1/2, which is the critical probability for edge percolation on the square lattice [5]. 
This motivates our next conjecture. 

Conjecture 5.1. For any e > 0, it holds that, almost surely, all components of R(TL*) are finite. 

Now we consider the question of whether percolation is a monotone property. We say that a 
point process LTo is transitive if, for every two nodes v, v' 6 LTo, there exists an isometry / : LTo — > LTo 
such that f(v) = v 1 . The open problem below concerns the question of whether transitivity is 
enough to obtain monotonicity in the percolation properties of balls moving as Brownian motion. 
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Figure 2: An example of a non-transitive configuration of balls that is not monotone. Each solid 
ball represents the superposition of 14 balls and white balls represent single balls. 

Question 5.2. Let Ilo be a transitive point process so that i?(IIo) is a connected set. Let lit be 
obtained from Ho by letting the nodes move as independent Brownian motions for time t. Then, 
if for some time to we have that R(Tlt ) has an infinite component almost surely, then, is it true 
that, for any t < to, R(Ht) also has an infinite component almost surely? Similarly, if for some t\ 
we have that it^II^ ) contains only finite components almost surely, then, does it hold that, for any 
t>t\, R(Tlt) also contains only finite components almost surely? 

Remark 5.3. We note that Question 0T2] above is false if we drop the condition that Ilo is transitive. 
For example, consider a tessellation of M? into squares of side length 6 and, in each square of the 
tessellation, consider the configuration of balls illustrated in Figure [21 where each ball has radius 
1/2, solid balls represent the superposition of 14 balls and white balls represent a single ball. It is 
easy to see that, at a sufficiently small time e, the union of the balls will not contain an infinite 
component almost surely. However, the density of balls is equal to 9xl 3 4 6 +18 = 4 and, as the balls 
move for a sufficiently large amount of time, their position will approach a Poisson point process 
which is known to percolate. 

5.2 Motion over graphs 

We now consider the case when the motion of the nodes is more restricted. First, let Ilo be the 
point process given by the integer points of R. For any node u £ Ilo, we let u + ( u (t) be its position 
at time t, where (Cu(t))t is a one-dimensional Brownian motion. Now, consider a sequence of m 
distinct nodes u±,U2, ■ ■ ■ ,u m such that B{ui, 1/2) and B(ui + \, 1/2) intersect for all i. We call such 
a sequence of nodes as a path. Let e be a sufficiently small positive constant, and consider only the 
nodes of Ilo whose displacement from time to time e is smaller than 1/2; we denote these nodes 
as good nodes. We claim that 

P (u\,U2, ■ ■ ■ , u m form a path at time e | Ui is good for all i) = — -. (15) 
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In order to see this, suppose, without loss of generality, that u\ < 112 < ■ ■ ■ < u m . For each node 
u € IIo, let Cu( e ) be the displacement of u from time to e given that u is a good node. Then, in 
order for B(ui+( Ul (e), 1/2) to intersect B(u 2 +Cu 2 ( e )i 1/2) we need that |^i+Cui( e ) — u<2 ~ C 2 ( e )l — ^ 
Since ui and U2 are good nodes, this condition translates to U2 + C 2 ( e ) — u i — d( e ) — 1j which 
in turn implies that (' U1 (e) > Cu 2 ( e )- Repeating this argument, we obtain the condition (' U1 (e) > 
C' U2 ( € ) — (u 3 ( € ) ^ ' ' ' ^ C m ( e )- Since the £' are independent and identically distributed, we have 
that P (C ui (e) > C 2 (e) > ™ • > C m ( e )) = V™ ! > which establishes jTgJ. 

We now consider a more general scenario. Let G be an infinite graph that is vertex transitive 
and has bounded degree. We assume that each edge of G has length 1, which gives a metric over 
G. Let IIo(G) be the point process given by putting one node at each vertex of G and define IL(G') 
as the point process obtained by letting the nodes of IIo(G) move for time t along the edges of G 
according to independent Brownian motions. Then R(Tlt(G)) is the union of balls centered at the 
nodes of lit and having radius 1/2 with respect to the metric induced by G. We note that the 
probability given in (fT5j) for any fixed path u\ , u% , . . . , u m of good nodes to form a path at a time 
e that is sufficiently small is at most 1/m!. This motivates our next conjecture. 

Conjecture 5.4. Let G be an infinite graph that is vertex transitive and has bounded degree. Then, 
for any t > 0, the region R(U t (G)) contains only finite components almost surely. 

5.3 Critical radius of point processes 

Here we let II be a point process over M? and consider the region R(H, r) as the union of balls of 
radius r centered at the nodes of IT. In this section, we only consider point processes with unit 
intensity and let ^(LT) be the smallest r for which R(Tl, r) contains an infinite component. It is 
intuitive to believe that point processes that are more organized have smaller critical radius; this 
is the core of our next conjecture. For more information on zeros of Gaussian analytic functions, 
we refer to [6]. 

Conjecture 5.5. Let IIl be any transitive point process with intensity 1 (as defined before Ques- 
tion \5.2\) . Let IIgaf be a point process given by the zeros of a Gaussian analytic function with 
intensity 1 and Hp be a Poisson point process with intensity 1. Then, 

r c (n L ) < r c (IIgaf) < r c (n P ). 

Finally, consider a Poisson point process II with intensity 1 over R d and let r c be the critical 
radius for percolation of balls centered at the nodes of II. Our last open problem concerns small 
perturbations of the critical radius. 

Question 5.6. Let e > and, for each node u G II, let X v be a uniform random variable over 
[—e, e]. For each node v G LT, add a ball of radius r c + X v centered at v. Will the union of the balls 
contain an infinite component almost surely? 
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A Standard large deviation results 

We use the following standard Chernoff bounds during our proofs. 

Lemma A.l (Chernoff bound for Poisson). Let P be a Poisson random variable with mean A. 
Then, for any < e < 1 , 

P(P> (l + e)A) <exp^-^(l-e/3)^ and P (P < (1 - e)A) < exp 

Lemma A. 2 (Chernoff bound for Binomial, see [TJ Lemma A. 1.4]). Let X\,X%, ■ ■ ■ ,X n be in- 
dependent Bernoulli random variable such that E [Aj] = pi. Let X = Y17=i-^i- Then, for any 
e > 0, 

P(X > (l + e)E[A]) < exp (- 2e2(E J X] H and P (P < (1 - e)A) < exp (-^- 



Lemma A. 3 (Gaussian tail bound [10^ Theorem 12.9]). Let X be a normal random variable with 



mean and variance a 2 . Then, for any R> a we have that P (X > R) < ,° exp ' 
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